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We study the 'spectral projector' method for the computation of the chiral condensate and the 
topological susceptibility, using iVV = 2+1 + 1 dynamical flavors of maximally twisted mass 
Wilson fermions. In particular, we perform a study of the quark mass dependence of the chiral 
condensate £ and topological susceptibility % top in the range 270 MeV < m K < 500 MeV and 
compare our data with analytical predictions. In addition, we compute % top in the quenched 
approximation where we match the lattice spacing to the Nf = 2 + 1 + 1 dynamical simulations. 
Using the Kaon, 77 and 77' meson masses computed on the Nf = 2 + 1 + 1 ensembles, we then 
perform a preliminary test of the Witten-Veneziano relation. 
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1. Introduction 

The Banks-Casher relation [ 1 ] connects the chiral condensate £ with the density of eigenmodes 
at the origin of the spectrum and thus with the infrared properties of the Dirac operator. The chiral 
condensate is obtained from the eigenvalue density p (X , in) in a triple limit, sending the volume V 
to infinity and the quark mass m as well as the eigenvalues X to zero (in this order), 

I 

— = lim lim lim p(X,m). (1.1) 

71 A->0m->0V->'°° 

Only when the eigenvalue density p(X,m) is non-zero at the origin, the chiral condensate 
does not vanish and hence the infrared properties of the Dirac operator are directly related to the 
mechanism of chiral symmetry breaking. 

One way to express p(X,m) is through the mode number v(M,ni), which is defined as the 
number of eigenmodes X of the considered Dirac operator squared below some cut-off mass M, 

v(M,m)=V f A dXp(X,m), A=y / M 2 -m 2 . (1.2) 

J-A 

The above considerations can -in principle- be taken over directly to the lattice as a way to 
compute the chiral condensate non-perturbatively. However, counting the eigenmodes below the 
cut-off M with M « O(100)MeV and taking the continuum limit more and more modes have to 
be taken into account for a fixed physical value of M. In fact, a direct counting of the low-lying 
eigenmodes is expected to show an 0{V 2 ) scaling behaviour and becomes prohibitively computer 
time expensive when the continuum limit is taken. 

Recently, however, a new method [2] to compute the mode number was developed, the so- 
called spectral projector method. The important advantage of this method is that it is computation- 
ally much faster than counting eigenmodes directly and scales only with the volume V. In addition, 
the concept of spectral projectors can be extended to evaluations of other quantities such as the 
topological susceptibility Xtop or the ratio of the pseudoscalar and scalar renormalization constants 
as explained in refs. [2, 3]. 

It goes beyond the scope of this proceedings contribution to detail the spectral projector 
method and we have to refer to refs. [2, 3] for a description of this method. The aim of this 
contribution is rather to see, how the spectral projector method works for computing the chiral 
condensate and the topological susceptibility in the case of the here used maximally twisted mass 
fermions formulation of lattice QCD [4]. In particular, in this work we are interested in the quark 
mass dependence of the chiral condensate and the topological susceptibility and we will work at 
only one value of the lattice spacing of a 0.0782 fm. All results are shown for a setup employing 
a mass-degenerate light quark doublet and a strange and a charm quark close to their physical val- 
ues, a situation we refer to as Nf = 2 + 1 + 1, see refs. [5, 6, 7] for simulation and analysis details. 
Employing several values of the quark mass will allow us to confront our data with predictions of 
chiral perturbation theory and to extract values for the chiral condensate in the chiral limit. We will 
also perform a first test of the Witten-Veneziano formula [8, 9] in this proceedings contribution by 
computing Kaon, tj and r\' masses on our dynamical Nj ■ = 2 + 1 + 1 configurations and the topo- 
logical susceptibility in the infinite quark mass limit (quenched approximation) matched, however, 
to the physical situation of the unquenched simulations. 
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2. Evaluation of the chiral condensate and topological susceptibility 

When computing the chiral condensate from spectral projectors, there are two important in- 
gredients. The first is a technical aspect: the spectral projectors are calculated from a stochastic 
estimate of the "inverse" of some suitable function of the lattice Dirac operator employed, see 
ref. [2]. Therefore, it needs to be investigated what is a sufficient number of stochastic noise vec- 
tors employed and how the stopping criterion for obtaining the solution of a Dirac equation needed 
to construct the spectral projector. The second aspect is more physical and originates from the fact 
that the chiral condensate is computed from the slope of the mode number as a function of the 
cut-off M. Before coming to our results, let us therefore briefly discuss the tests we have made for 
both issues. In the following, we will use M* as the cut-off for the mode number counting. M* is of 
a very similar size as M and plays the role of an adjustable parameter to optimize the simulations. 

As a very first test, we performed a comparison between the explicitly computed mode number 
and the values from the spectral projectors. In this test we found a perfect agreement demonstrating 
that our implementation of the spectral projectors in the tmLQCD package [10] is correct. 

Looking at the mode number itself as function of M, see fig. 1 (left), we indeed can identify 
a linear behaviour of the mode number which will allow us in the following to extract the chiral 
condensate. We also show in fig. 1 (right) our results for changing the number of stochastic sources 
and the (relative) stopping criterion. As a conclusion from this study we found that with already 6 
stochastic sources the corresponding error saturates, nevertheless we used 8 sources to remain safe. 
In addition, we found that the stopping criterion can be chosen rather loosely and even a choice 
of 10~ 2 gave completely consistent results. Nevertheless, for our work we decided to choose a 
stopping criterion of 10~ 6 to be on the safe side. 
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Figure 1: (left) The mode number as a function of M - from explicit computation of eigen- 
modes (line) and from spectral projectors (points), (right) The influence of the number of stochastic 
sources and relative precision of solving the Dirac equation. 

2.1 Results 

After the tests described in the previous section, we proceeded to compute the mode number 
from spectral projectors. As said above, for this work we used only one value of the lattice spacing 
of a « 0.0782fm corresponding to j8 = 1.95, as determined by ETMC [6]. While working here 
only at one value of the lattice spacing, we have computed the mode number at several values of the 
quark mass, see tab. 1, where we give the bare light twisted mass parameter in lattice units, as well 
as the pion mass. Note that for all the parameters shown in tab. 1 the theory was tuned to maximal 
twist. Our typical statistics for computing the mode number and the topological susceptibility has 
been 200 configurations that were separated by 20 HMC trajectories. 
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lattice 




M n (MeV) 


1.95 


32 3 x 64 


0.0025 


270 


1.95 


32 3 x 64 


0.0035 


320 


1.95 


32 3 x 64 


0.0055 


390 


1.95 


32 3 x 64 


0.0075 


455 


1.95 


24 3 x 48 


0.0085 


490 



Table 1 : Parameters of ensembles used to computed 
the mode number and the topological susceptibility. 
We give the light bare twisted mass parameter a\i\ in 
lattice units as well as the approximate pion masses. 
The value of j8 = 1.95 yields a a lattice spacing of 
a «0.0782fm, seeref. [6]. 



2.2 Evaluation of the chiral condensate 

We have computed the chiral condensate £ through a study of the average number of eigen- 
modes of D'D with X < (M*) 2 , computed using the spectral projector method. 

Once the mode number is computed at several values of M*, the chiral condensate can be 
calculated directly from the derivative of the mode number with respect to M* [2], 




Mi 



dM 



-v R (M*,n, iR ). 



(2.1) 



R / R 

In eq. (2.1), fi^p and denote the renormalized twisted mass and cut-off parameters. The 
renormalization constant Zp (MS, \l=2 GeV) needed for obtaining these renormalized quantities 
has been computed by ETMC in a dedicated four flavour simulation, see ref. [11]. Since the mode 
number is renormalization group invariant (i.e. Vp{Mp,pLip) = v(M,/i/)) [2], from eq. (2.1) we 
hence obtain directly the renormalized chiral condensate at a scale that is inherited from Zp. 
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Figure 2: (left) The mode number as a function of Mp\ for ~ 14 MeV and the corresponding 
linear fit. (right) The chiral condensate as a function of the renormalized quark mass. The 
straight line indicates a linear extrapolation of Lp to the chiral limit. 



In fig. 2 (left) we show an example of a behaviour of the mode number as a function of M^, for 
a renormalized quark mass of 14MeV. For the four values of we have used (and similarly as in 
fig. 1), we observe a linear behaviour of the mode number in the range 60MeV < < 120MeV. 
Such linear behaviour of v in a range of comparable values of Mp\ was also observed in ref. [2]. 
From this linear behaviour, which we see at all five quark masses employed, we can extract the 
renormalized chiral condensate Lp using eq. (2.1) at a given value of the renormalized M*. In 
fig. 2 (right) we show as a result the depence of Lp on the renormalized quark mass. Extrapolating 
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L R linearly to the chiral limit, we rind L l R /3 (MS,n = 2GeV) = 312(l)(13)MeV. The first error 
is purely statistical and the second originates from the uncertainty of the renormalization constant 
Z P (MS,n = 2GeV) = 0.462(13) [12]. 

2.3 Topological Susceptibility 

As another quantity accessible to the method of the spectral projectors, we have computed 
the topological susceptibility % top following ref. [3]. In this reference it was demonstrated that 
besides the bare % top also the renormalized one can be obtained by solely using observables defined 

through spectral projectors. Since we know the necessary renormalization factor, the ratio of the 

z 2 

scalar to the pseudoscalar renormalization constants ^ available to us from the Nf = 4 simulations 
of ETMC [12], we decided to only compute the bare value of Xtop from spectral projectors and to 
perform the renormalization using the results from ref. [12], i.e. ^ = 0.685. 

In this way, we have computed % top at five values of the bare quark mass, listed in tab. 1. 
Before showing our results, we remark that we have performed a test of the dependence of Xtop as 
a function of the cut-off parameter M R . We found that in the range 90MeV < M R < 130MeV Xtop 
is constant as a function of M R and we therefore decided to use a value of M R « lOOMeV for the 
computation of all subsequent values of % top . 
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Figure 3: The chiral behaviour of the topolog- 
ical susceptibility with 5 values of the quark 
mass. The linear fit represents the tree-level for- 
mula of #PT, Xtop = and yields a value 
of the chiral condensate. 
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Figure 4: The critical hopping parameter K^t 
for the quenched ensemble at j8 = 2.67 as a 
function of fii. The linear fit yields the chiral 
limit value of fQ r it which is used to realize max- 
imal twist. 



Fig. 3 shows the chiral behaviour of the topological susceptibility. We have fitted our re- 

1 /3 

suits with the tree-level formula of chiral perturbation theory, Xtop = HiE/2. The fit yields L R = 
282(5) (13) MeV, where again the systematic error is dominated by the uncertainty in the renor- 
malization constants. 

3. Witten-Veneziano formula 



As said above, with spectral projectors we have a method at hand that allows for a rather cheap 
computation of the topological susceptibility, which is, in addition, well defined, i.e. it does not 
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suffer from short distance singularities. Moreover, twisted mass fermions are advantageous for 
computing disconnected (singlet or OZI) quantities as, e.g., the masses of the v\ and v\' mesons 
(see the contribution of V. Drach to this conference and refs. [13, 14]). 

Thus, it is very tempting to attempt a non-perturbative test of the Witten [8] - Veneziano [9] 
formula, which provides an elegant explanation for the origin of the unexpectedly large mass of the 
tj' meson. The Witten- Veneziano (WV) formula relates the masses of the Kaon, tj and tj' mesons 
to the topological susceptibility x°° where the oo index reminds us that the topological susceptibility 
needs to be computed at infinite quark mass, i.e. in the pure gauge theory. The formula then reads, 

f 2 

±^( m \+m\,-2m\)=l^ (3.1) 

where f n is the pion decay constant. Eq. (3.1) cannot be derived in a rigorous way, but some (mild) 
assumptions are required. It is therefore most worthwhile to test the formula by direct and non- 
perturbative lattice simulations, since the WV formula is one of the most fundamental relations in 
QCD and clearly points out the importance of topology. 



3.1 Strategy to compute %„ 

In order to test the WV formula (3.1), the topological susceptibility x«> needs to be computed 
in the pure gauge theory, but at a matched physical situation. To fulfil this condition, we performed 
a scan in /3 using the same (Iwasaki) gauge action as used in [6] and searched for the value of j8 that 
gives the same value of the force parameter ro as the (chirally extrapolated) one of our dynamical 
ensembles at j3 = 1.95. As a result, we obtained jS = 2.67. The meson masses that are needed in 
the WV formula were evaluated for a bare twisted mass parameter of \L\ = 0.0055. 

Since the topological susceptibility from spectral projectors is a fermionic quantity and we 
want to use maximally twisted mass fermions, a first step is the tuning to maximal twist. This 
amounts to tuning the bare Wilson quark mass, or equivalently the hopping parameter K, to its 
critical value, m cr - lt , leading to K^t = 1/(8 + 2m cr j t ). Following the strategy introduced in ref. [15], 
we have computed K^ r ;t at different values of the twisted mass /I/ by demanding that the PCAC 
quark mass vanishes. In the end, we have performed a chiral extrapolation letting /I/ approach 
zero. Our critical value of K^t is then the one in the chiral limit. The results of this procedure are 
shown in fig. 4, which shows that a linear extrapolation to zero twisted mass parameter is justified. 

Since in the pure gauge theory and with the Iwasaki gauge action the renormalization constants 
Z$ and Zp are not available to us, we decided to follow in this case refs. [2, 3] and compute the 
renormalized x°° solely from suitable expectation values employing spectral projectors. 

In this way, we finally obtained a value of x°° which is listed in tab. 2 together with our 
results of the meson masses and f K relevant for the WV formula. Putting everything together and 
multiplying both sides of eq. (3.1) with to make it dimensionless, we find for the left hand side 
of the WV formula 0.036(8) and for the right hand side 0.053(18). Although within the errors the 
WV formula is fulfilled, our present data clearly do not allow us to perform a stringent test. 



am-q 




aniK 


afn 


« 4 Z~ 


0.230(10) 


0.384(24) 


0.2280(4) 


0.0656(2) 


0.000050(17) 


r 4 oS} { m2 r i + m n' ~ lm "k) = °- 036 ( 8 ) 


t$Xoo = 0.053(18) 



Table 2: Results of the 
meson masses and in 
pure gauge theory. 
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4. Conclusion 

In this proceedings contribution we have explored the potential of the newly introduced spec- 
tral projectors method when applied to maximally twisted mass fermions in the Nf = 2 + 1 + 1 
setup. Our analysis has used only one value of the lattice spacing, but several quark masses which 
allowed us to compute the chiral condensate in the chiral limit, both from the quark mass depen- 
dence of the condensate itself and the topological susceptibility. In addition, we can compare these 
values to the ones of ref. [6], where the chiral condensate has been extracted from chiral perturba- 
tion theory fits to the pion mass and decay constant. In tab. 3 we show these different results. The 
agreement between the extraction of £ using very different methods is reassuring. 



4 /3 (M =2GeV)[MeV] 


spectral proj 


chiral fits 


Xtop 


312(1)(13) 


276(2)(11) 


282(5)(13) 



Table 3: Results for the chiral condensate ob- 
tained with 3 different methods. 



As another step, we have performed a first test of the Witten-Veneziano formula. Unfortu- 
nately, at the moment our accuracy does not allow for a stringent test of this fundamental relation 
between meson masses and the (quenched) topological susceptibility. However, the fact that we 
found reasonable errors already is quite promising that in the future a more precise test can be 
performed. Clearly, in this work a number of systematic effects could no be considered yet. In 
particular, it will be interesting to understand the effects of the lattice spacings. 

As a last remark, we want to mention that the spectral projector expectation values used to 
compute the topological susceptibility show a high sensitivity of autocorrelations stemming from 
topology. Since these expectation values are rather cheap to compute, it can thus be envisaged that 
such quantities can be used to scrutinize simulations of lattice QCD for possible large autocorrela- 
tion times. 
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